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Abstract. This article considers a postertori error estimation of specified function- 
als for first-order systems of conservation laws discretized using the discontinuous 
Galerkin (DG) finite element method. Using duality techniques, we derive exact 
error representation formulas for both linear and nonlinear functionals given an as- 
sociated bilinear or nonlinear variational form. Weighted residual approximations of 
the exact error representation formula are then proposed and numerically evaluated 
for Ringleb flow, an exact solution of the 2-D Euler equations. 


1 Introduction 

A frequent objective in numerically solving partial differential equations is 
the subsequent calculation of certain derived quantities of particular inter- 
est. e.g., aerodynamic lift and drag coefficients, stress intensity factors, etc. 
Consequently, there is considerable interest in constructing a posteriori error 
estimates for such derived quantities so as to improve the reliability and ef- 
ficiency of numerical computations. For an introduction to o posteriori error 
analysis see Eriksson et al. [9], related work by Estep et al. [13], Parashivoiu 
et al. [15], and the recent report of Oden and Prudhomme [14]. For hyperbolic 
problems and applications in fluid mechanics see Johnson et al. [12], Giles et 
al. [10], Becker and Rannacher [4] and Siili [16]. 

This article revisits the topic of a posteriori error estimation of prescribed 
functionals with special emphasis and consideration given to nonlinear sys- 
tems of conservation laws discretized using the discontinuous Galerkin (DG) 
finite method, see for example Johnson and Pitkaranta [11], Bey and Oden 
[5], and Cockburn et al. [7,8]. In a departure from this previous work, our 
DG formulation for systems of conservation laws uses entropy symmetriza- 
tion variables as discussed in detail in the companion papers by the second 
author [3.2.1]. 

In Section 2, we briefly Feview the abstract model for a posteriori error 
estimation of functionals. Next, we consider the DG method for nonlinear 
systems of conservation laws and derive concrete error estimates in terms 
of element residual and weight formulas. Section 4 numerically assesses the 
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sharpness of these estimates for the specific example of Ringleb flow which 
has a known exact solution via hodograph transformation. 


2 .4 Posteriori Error Estimation of Functionals 


Abstract model problem. In this section, we give an abstract presentation 
of a posteriori error estimation for functionals based on duality techniques. 
Consider the following abstract variational problem: find u € X such that 

A(g;u,v) = 0 VueX (2.1) 

and the corresponding discrete problem: find u k € X h such that 

A(g; u h , v h ) = 0 Vu A e X h . (2.2) 

Here X is a suitable function space, X h C X is a discrete space, for instance, 
discontinuous piecewise polynomials of degree k, and g some prescribed data. 
Note that boundary conditions are weakly imposed in the variational state- 
ment thus permitting both u and v to reside in X. For brevity, we sometimes 
write A(uh. v /,) = A(g:Uh,Vh)- Our objective is to estimate the error 

M{u) - M(u h ), (2.3) 

in a given functional M (*). The first step is to derive an error representation 
formula. 


Error representation : linear case. We first assume that A(-, ■) and M( ) are 
both linear. To derive a representation formula for the error (2.3), we intro- 
duce the dual problem: find such that 


A(v, <P) = M(v) V v € X. 
Setting v = u - u h in (2.4) yields 


(2.4) 


*V(tt)-.\/(li fc ) = -J/(w-tl fc ) 

= A(u - u h ,$) 

= A(U - - 7T*#) 

~ A(u,<P - 7T/,#) - A(Uh,$ - 7T/,#) 

= -A(u h ,$ - n h $) 


(linearity of M) 
(2.4) 

(orthogonality) 
(linearity of A) 
( 2 . 1 ), 


where £ Xh C A is an interpolant of 4L Thus we have the error repre- 
sentation formula 


.\/(u) - M{u h ) = - A{g,Uh,<P - 7 r^#). 


(2.5) 
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3 A Posteriori Error Estimates for the DG Method 

First-order nonlinear system. Consider the prototype conservation law prob- 
lem: find u : fl -> IR m such that 


L{u) = =0 in f?. (3.1) 

A~(n; u)(g — n) — 0 on F. 

where Q is a domain in lR d with boundary r with local exterior normal vector 
n and A{n: u) = n,/‘ u is the flux Jacobian matrix. In addition, A(n;g,u) 
denotes the mean value matrix obtained from the path integration 

.4(n;M)= I' A(n;t+G(s-t)) dO (3.2) 

Jo 

and P ± {n:s,t) the associated characteristic projectors. Throughout, we as- 
sume that u denotes the symmetrization variables so that the matrices .4 are 
necessarily symmetric. 

Next, consider a finite element tessellation T of Q composed of nonover- 
lapping elements T t . T = UT; t T t fl Tj = 0, i / j and Ft the tessellated 
boundary. The prototype system can be restated in variational form 1 * : find 
u € A' such that V v G X 

A(g : u.v) = ^2 ((£(u), v)t 4- (*4“(n;g,u)(5 - u),v)drnr T ( 3 * 3 ) 

T €T 

+ ( P~(n ; u-, «+)[/(«; u)]t, V-)3T\r r ) 

Note that other mathematically equivalent formulations are possible by group- 
ing together terms element-wise and edge-wise. The above particular grouping 
has been chosen as it reflects a discrete balance of conserved quantities on an 
element-by-element basis. In Section 4, we briefly revisit the possibility of al- 
ternate groupings although our numerical results show' that the element-wise 
grouping presented above yields superior estimates. 


A posteriori error estimate residuals and weights. A straightforward appli- 
cation of Cauchy-Schwarz inequality (with Ao introduced from entropy sym- 
metrization theory for dimensional consistency) (3.3) yields the following 


1 In actual implementations it is desirable to use an integrated-by-parts form (see 

for example [3,2]) so that exact discrete conservation is achieved on elements with 

inexact quadrature and/or inexact path integration (3.2). 
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Error representation, nonlinear case. Consider now the case of a nonlinear 
variational form -4(-. •) and functional M(-). To perform the analysis given 
above, we introduce the following mean value linearizations 

,4(g; u. v ■ = A(g: u h . v) + A(g. u h , u; u - u h . v) V v € X (2.6) 
M(u = M(uh) + M(Uh, u; u — u/,). (2.7) 

and the dual linearized problem: find <P e X such that 

■Alg, Uft,u:w,#) = M(uh,u;v) V » 6 I. (2.8) 

In addition, we have the following orthogonality relation 

A(g,u h ,u;u - u h ,v h ) = 0 V v h € X h . (2.9) 

Proceeding in the same fashion as above, using simplified notation for brevity, 


.\/(ti) — M(uh) = A/(u — u/,) (2.7) 

= A(u - u h ,$) (2.8) 

= A(u - u h , # - tt/,#) (2.9) 

= A(u,<P - n h <t>) - A(u h ,$ - ir h $) (2.6) 

= -A{u h .$ - 7r h #), (2.1) 

thus yielding the following final error representation formula 

M{u) - \I{u h ) = -A{g;u h ,$ - w h $). (2.10) 


Abstract a posteriori error estimates. Starting from (2.5) or (2.10), we derive 
various error estimates by estimating the right hand side of (2.10) using 
standard inequalities. Later, the sharpness of these inequalities is numerically 
assessed. Consider the following sequence of direct estimates 

|A/(u) - A/(u h )| = 1 (</;«*,<? - 7r*<?)| (2.11) 

T 1 

^ 5Zl-4T(ff;u h ,#-jr A #)| (2.12) 

T 

^^2 R r{u h ) (2.13) 

T 

where Aj{-.) denotes the restriction of A(-,-) to the element T. Further 
ctr(w/.) is a computable estimate of the residual of u h on T. and W T (<P) is 
a weight on T describing the local influence of #, both are defined below. 
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o 


element residuals Rt and weights W T for use in (2.13) 


R T {u h ) 


H £ («*)lb 0 -.r 

P (n: ^)[/(n; !/*)]_ || QT\r r 


M {n.g:u h ){g -u h )\\ x -i rT 




' ll*-^*IU 0i r ' 

11 ^ - dT\r r 

L Ut-K^hoTT \ 


(3.4) 


(3.5) 


Approximating the dual pro6/em. The weight formulas (3.5) require the cal- 
culation of the quantity # — 7 from the dual problem w T hich requires a 
priori knowledge of both u and u ti for use in the mean value linearizations 
(2.6) and (2.7). Since u is not generally known, we supplant this calculation 
with the approximate discrete counterpart to (2.8): find <Ph € X * such that 

A{g,u h ,u k \v h ,$ h ) = Yl{n h) u h \v h ) V v e X h . (3.6) 

Observe that A(g, u^; v, <Ph) and M{uh, u^; Vh) are precisely the Jaco- 
bian linearized forms of the respective operators. Using the techniques de- 
scribed in Barth [3.2], exact Jacobian derivatives of the DG scheme for sys- 
tems of conservation laws have been derived and used in all subsequent cal- 
culations. We have investigated the computation of the needed dual solution 
terms using two different techniques: 

(1) High-order approximation. Suppose denotes a numerical solution 

computed in Embed u in l > k and approximate # — 71 ^# % 

k)< &h ) - This technique is employed in the calculations given below. 

(2) Recovery post-processing. Let : x[ k) ^ xff } denote a recovery 

operator, approximate l > k. Recovery operators 

based on local compact supported least-squares fitting are considered in a 
forthcoming report by the present authors. 


4 Numerical Results 

To evaluate the accuracy of the error representation formulas given in Sect. 2, 
Ringleb flow (an exact solution of the 2-D Euler equations obtained via hodo- 
graph transformation, see [6]) is computed in the channel geometry shown in 
Fig. 4.1(a). Next, the vertical force component exerted on the channel walls 
is computed from the functional 

Mxp{u) = f {'&n)p{u)dl 
Jr w^n 


(4.7) 
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fa) Coarsest tnangulation 
(342 elements). 



(b) Primal solution iso- 
density contours. 



(c) Dual solution iso- 
density contours. 


Fig 4.1. Ringleb dow test problem. Pr.mal and dual solutions calculated using the 
DG discretization with cubic elements for the vertical force funct.onal A/,(u). 


with pin) the fluid pressure and 9 a constant vertical vector. k<>density 
contours of the Ringleb primal and dual solutions are given in lgs. - b 
We now evaluate the validity and sharpness of the error es.tma.e formulas 
(9 ii)-(2.13) and (3.5). In Fig. 4.2 we graph for constant (a) and linea 
i approximation: (o| the exact error; <x) estimate (2.11), W — 
(2 12)' (7) estimate obtained from element-edge form of (3.3), (□) estuna 
,3.5) In all case, the dual problem is defined by (3.6) and solved using cubic 
polynomials. The difference between (o) and (x) is caused by lineanzatio 
ce replacing « in (2.8) by u„ to get (3.6). This appears to be a very snud 
error. Next, the more significant loss due to use of the Tnang e nequ l y i 
crraphed in (A). This prevents cancellation between elements. Further error 
fs introduced via Cauchy-Schwarz (□) thus preventing cancellation within 
the element. Finally, note that the element based estimate (A) is ^ 
superior to the element-edge based estimate (V), where in the latter case 
contributions are grouped together in such a way that dement consermt^n 
is violated. Based on our numerical experimentation, we propose the adaptive 

method: ... 

• Evaluate a stopping criterion via |.4(u h ,4‘ - * h9 ' n '' 

• Evaluate an adaptation criterion via |x4(u/,, 9 - vta<P)t|- 

In addition, the adaptation criterion may be further improved by the use 
of sharpened variants of the Triangle and Cauchy-Schwarz Inequalities. We 
consider these topics further in a forthcoming paper. 
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